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Abstract 

We search for lepton- flavor- violating r decays into a lepton (electron or muon) and an /o(980) 
meson using 671 fb" 1 of data collected with the Belle detector at the KEKB asymmetric-energy 
e + e~ collider. No events are observed and we set the following 90% C.L. upper limits on the 
branching fraction products: B{t~ -> e - / (980)) x £(/ (980) -> Tr+Tr") < 3.2 x 10~ 8 and B(t~ -» 
A i_ /o(980)) x £>(/o(980) — ► ir + ir~) < 3.4 x 10~ 8 . This is the first search performed for these modes. 

PACS numbers: 11.30.Fs; 13.35.Dx; 14.60.Fg 
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INTRODUCTION 



Lepton flavor violation (LFV) in charged lepton decays is forbidden in the Standard Model 
and remains highly suppressed even if neutrino mixing is included. However, LFV is expected 
in various extensions of the Standard Model, such as supersymmetry, leptoquark and many 
other models, e^g., various MSSM scenarios consider r~ lepton decays into semileptonic 
final states [l|, |2|. If LFV occurs via a Higgs mediated LFV mechanism, r leptons can 
decay into £~/o(980), where i~ is either an electron or a muon, through a scalar Higgs 
boson 0[t]. The decays £~ir°, i~r\ and £~r( are mediated by a pseudoscalar Higgs boson 
while can be mediated through both scalar and pseudoscalar Higgs bosons 

Searches for r~ — > £~/ (980) have not yet been reported while there are already upper 
limits on the branching fractions for r _ — > £~ir (77, 7/) and r~ — > £~fi + fi~ at the O(10~ 8 ) 
level (5), 0] . With certain combinations of new physics parameters the branching fractions for 
r~ —¥ £~/ (980) decays can be as high as 1CT 7 , which is already accessible in high- statistics 
i?-factory experiments. 

Here we report our search for r~ — > i~ / (980) with a data sample of 671 fb _1 collected at 
the T(4S) resonance and 60 MeV below it with the Belle detector at the KEKB asymmetric- 
energy e + e~ collider [zj]. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals, all located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside the coil is instrumented to detect K° L mesons and to identify muons 
(KLM). The detector is described in detail elsewhere j8|. 

We use particle identification likelihood variables based on the ratio of the energy de- 
posited in the ECL to the momentum measured in the SVD and CDC, the shower shape 
in the ECL, the particle range in the KLM, the hit information from the ACC, the dE/dx 
information in the CDC, and the particle time-of-flight from the TOF. To distinguish hadron 
species, we use likelihood ratios, for instance, V(i/j) = £i/(£% + Cj), where Ci (Cj) is the 
likelihood for the detector response to the track with flavor hypothesis % (J). For lepton 
identification, we form likelihood ratios Vie) jsj and V(fi) jlo| based on the electron and 
muon probabilities, respectively, which are determined by the responses of the appropriate 
subdetectors. 

In order to estimate the signal efficiency and to optimize the event selection, we use 
Monte Carlo (MC) samples. The signal and the background events from generic t + t~ 
decays are generated by KKMC/TAUOLA [ll|. In the signal MC, we generate r + r~, where 
a t decays into a lepton and an /o(980) meson according to two-body phase space, and the 
other r decays generically. In the signal MC the /o(980) meson decays into 7r + 7r~ with mass 



and width parameters of the /o(980) based on the result of our study of 77 — > tx + ti~ [12 
Other backgrounds, including BB and e + e~ — > qq (q = u,d,s,c), Bhabha and two-photon 
processes are generated by EvtGen BHLUMI [3], and AAFH 15], respectively. All 



kinematic variables are calculated in the laboratory frame unless otherwise specified. In 



[|] Throughout this paper, charge-conjugate modes are implied unless stated otherwise. 
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particular, variables calculated in the e + e center-of-mass (CM) system are indicated by 
the superscript "CM". 



EVENT SELECTION 

We search for t + t~ events in which one r (the signal r) decays into a charged lepton 
and two charged pions from an /q(980) meson decay and the other r (the tag r) decays into 
a final state with one charged track and any number of additional photons and neutrinos. 
Candidate r-pair events are required to have four tracks with zero net charge. 

We reconstruct four charged tracks and any number of photons within the fiducial volume 
defined by —0.866 < cos0 < 0.956, where 9 is the polar angle relative to the direction 
opposite to that of the incident e + beam in the laboratory frame. The transverse momentum 
(pt) of each charged track and the energy of each photon (E 7 ) are required to satisfy p t > 
0.1 GeV/c and E 1 > 0.1 GeV, respectively. For each charged track, the distance of the 
closest point with respect to the interaction point is required to be less than 0.5 cm in the 
transverse direction and less than 3.0 cm in the longitudinal direction. 



Using the plane perpendicular to the CM thrust axis [la ], which is calculated from the 
observed tracks and photon candidates, we separate the particles in an event into two hemi- 
spheres. These are referred to as the signal and tag sides. The tag side contains one charged 
track while the signal side contains three charged tracks. We require the charged track on 
the signal side to be identified as a lepton. The lepton identification criteria are V{£) > 0.95 
and momentum greater than 0.8 GeV/c. The electron (muon) identification efficiency is 
92% (86%) while the probability to misidentify a pion as an electron (a muon) is below 0.2% 
(2%). In order to take into account the emission of bremsstrahlung photons from electrons, 
the momenta of photons within 0.05 rad of the electron tracks are added to the momentum 
of the electron candidate track. To reduce generic r + r~ and qq background events, we veto 
events with a photon in the signal side. 

To ensure that the missing particles are neutrinos rather than photons or charged particles 
that pass outside the detector acceptance, we impose requirements on the missing momentum 
Pmiss; which is calculated by subtracting the vector sum of the momenta of all tracks and 
photons from the sum of the e + and e~ beam momenta. We require that the magnitude 
of Pmiss be greater than 0.8 GeV/c, and that its direction point into the fiducial volume of 
the detector. Furthermore, we reject the event if the direction of the missing momentum 
traverses the gap between the barrel and the ECL endcap since an undetected photon may 
fake a neutrino candidate(s). 

The / (980) meson is reconstructed from two oppositely-charged tracks assuming a pion 
mass for both tracks in the signal side. The charged pion from the /o(980) decay is identified 
by a condition V(K/n) < 0.6 for two tracks except a lepton candidate in the signal side. 
Figure [T] shows the invariant mass of 7r + 7r~ in the signal side. We require that the invariant 
mass satisfy the condition 906 MeV/c 2 < M n + n - < 1065 MeV/c 2 . In order to avoid fake 
/o(980) candidates from photon conversions (i.e. 7 — > e + e~), we require that two pion 
candidates have V(e) < 0.1. Furthermore, we require that V(p) < 0.1 to suppress two- 
photon e + e~/i + /i~ background events. We require that the cosine of the opening angle 
between the lepton and the reconstructed / (980) candidate in the CM system {cos8^^ aa _f ) 
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FIG. 1: I The invariant mass distribution of ir + ir~ in the signal region. While the signal MC (r~ — > 
^~/o(980)) distribution is normalized arbitrarily, the data and background MC are normalized to 
the same luminosity. Selected regions are indicated by the arrows from the marked cut boundaries. 



be 0.7 < cos ^ on _ /o < 0.98. 

To reject qq background, we require that the magnitude of thrust (T) be in the range 
0.90 < T < 0.97 (see Fig. i (a)). We also require 5.29 GeV < < 10.00 GeV, where 
E™ is the total visible energy in the CM system, denned as the sum of the energies of the 
lepton, / (980) candidate, the charged track in the tag side (with a pion mass hypothesis) 
and all photon candidates (see Fig. [2] (b)) . 

To suppress the BB and qq background, it is required that the number of photons in the 
tag side n^ AG be n^ AG < 3 if the charged track is the tag side is an electron or a muon (a 
leptonic tag) and n^ AG < 2 if the track in the tag side is a hadron (a hadronic tag). 

Since neutrinos are emitted only in the tag side, the direction of p miss should be on the 
tag side of the event. The cosine of the opening angle between p m i SS and the charged track in 
the tag side in the CM system, cos0^_ miss , should be in the range 0.4 < cos0£^_ miss < 0.98 
(see Fig. [2] (c)). The reconstructed mass of a charged track (with a pion mass hypothesis) 
and photons in the tag side m tag , is required to be less than 1.00 GeV/c 2 (see Fig. [2] (d)). 
For all kinematic distributions shown in Fig. [21 reasonable agreement between the data and 
background MC is observed. 

Finally, to suppress backgrounds from generic t + t~ and qq events, we apply a selection 
based on the magnitude of the missing momentum p m i SS and the missing mass squared m^^. 
We apply different selection criteria since the number of emitted neutrinos is two for the 
leptonic tag while it is one for the hadronic tag. We require the following relation between 
p miss and m 2 miss : p miss > -5.0 x m 2 miss and p miss > 7.0 x m 2 miss - 1.0 for the hadronic tag and 
Pmiss > —7.0 x m^ iss — 1.0 and p miss > 1.8 x m^ iss — 0.8 for the leptonic tag, where p m i SS is in 
GeV/c and m miss is in GeV/c 2 (see Fig. [3]). While this cut retains 56% of the signal, 97% 
of the generic t + t~ and 96% of qq backgrounds are removed. 
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FIG. 2: Kinematic distributions used in the event selection: (a) the magnitude of the thrust vector 
using all track and photon candidates; (b) total energy using all track and photon candidates in 
the CM system; (c) the cosine of the opening angle between a charged track in the tag side and 
missing particles in the CM system (cos 0^_ miss ); (d) the reconstructed mass in the tag side using 
a charged track and photons. While the signal MC (t~ — > /z~/o(980)) distribution is normalized 
arbitrarily, the data and background MC are normalized to the same luminosity. Selected regions 
are indicated by the arrows from the marked cut boundaries. 



SIGNAL AND BACKGROUND ESTIMATION 

The signal candidates are examined in the two-dimensional plot of the l~ /o(980) invari- 
ant mass (Mef ) versus the difference of their energy from the beam energy in the CM 
system {AE). A signal event should have M^ close to the r-lepton mass (m T ) and AE 
close to zero. For all modes, the M^ and AE resolutions are parameterized from fits 
to the signal MC distributions, with an asymmetric Gaussian function that takes into ac- 
count initial-state radiation. These Gaussians have widths cr^l^ low = 4.8/5.6 MeV/c 2 and 
^high/low = 12 ^ ig MeV for the r _ _^ e -/ (980) mode, and a^ gh / low = 5.1/5.3 MeV/c 2 and 



a AE = 13/19 MeV for the r — > ji / (980) mode, where the "high/low" superscript 
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FIG. 3: Scatter-plots of p m i ss vs. m ^ m : (a) and (b) show the signal MC (t _ — ► ^~/ (980)) and 
the generic t + t~ MC distributions, respectively, for the hadronic tags while (c) and (d) show the 
same distributions for the leptonic tags. Selected regions are indicated by lines. 



indicates the higher/lower side of the peak. 

To evaluate the branching fractions, we use elliptical signal regions that contain 90% of 
the MC signal events satisfying all selection criteria. We blind the data in the signal region 
until all selection criteria are finalized so as not to bias our choice of selection criteria. 
Figure H] shows scatter-plots for the data and the signal MC distributed over ±20cr in the 
Mif —AE plane. For the e~/o(980) mode, the dominant background comes from two-photon 
processes in the AE > GeV region, while qq and generic t + t~ events are negligible due 
to the low electron fake rate. For the /x - /o(980) mode, the dominant sources of background 
are random combinations of three pions from qq or generic t + t~ processes in which one 
pion is misidentified as a muon. The qq events are distributed in the AE > GeV region 
while t + t~ events are distributed in the AE < GeV and M M j < m T region. We estimate 
the expected number of background events in the signal region using the MC distribution 
with looser selection criteria and extrapolation from the data in the sideband region denned 
as the ±20cr region without the signal ellipse. The signal efficiency and the number of 
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TABLE I: The signal efficiency^), the number of expected background events (-/Vbg) estimated 
from the sideband data, total systematic uncertainty (<7 sy st), the number of observed events in the 
signal region (iV bs), 90% C.L. upper limit on the number of signal events including systematic 
uncertainties (sgo) and 90% C.L. upper limit on the branching fraction including the branching 
fraction /q(980) — > tt + tt~ (JAC) for each individual mode. 



Mode 


e{%) 




O-gyst (%) 


Aobs 


«90 


^£(10" 8 ) 


T~ -> e-/ (980) 


5.80 


0.10 ±0.07 


11.5 





2.41 


3.4 


t- -> ^~/o(980) 


6.02 


0.11 ±0.08 


10.8 





2.40 


3.2 



expected background events for each mode are summarized in Table [B After estimating the 
background, we open the blinded regions and find no candidate events in either mode. 

The dominant systematic uncertainties for this analysis come from / (980) selection and 
tracking efficiencies. The systematic uncertainty for / (980) selection is studied by changing 
mass and width of the / (980) using uncertainties from the result of our study 12j, and 
is 9.9% and 9.2% for e~/ (980) and /i _ / (980), respectively. The uncertainty due to the 
charged track reconstruction is estimated to be 1.0% per charged track, therefore, the total 
uncertainty due to the charged track finding is 4.0% for both modes. The uncertainties due 
to lepton identification are 2.2% and 1.9% for electrons and muons, respectively. The uncer- 
tainty due to pion identification is 1.0% per pion. The uncertainties due to MC statistics 
and luminosity are estimated to be (2.7-2.9)% and 1.4%, respectively. The uncertainty due 
to the trigger efficiency is negligible compared with the other uncertainties. All these uncer- 
tainties are added in quadrature, and the total systematic uncertainties for the e~/ (980) 
and /i _ / (980) mode are 11.5% and 10.8%, respectively. 



UPPER LIMITS ON THE BRANCHING FRACTIONS 

We set upper limits on the branching fractions of r~ — > £~f Q (980) based on the Feldman- 



Cousins method 17j. The 90% C.L. upper limit on the number of signal events including a 



systematic uncertainty (sgo) is obtained using the POLE program without conditioning 
based on the number of expected background events, the number of observed events and the 
systematic uncertainty. The upper limit on the branching fraction (B) is then given by 

B(t~ -> r/„(980)) x fl(/ (980) - tt+tT) < (i) 

where N TT is the number of T + r~pairs, and e is the signal efficiency. The value N TT = 616.6 x 
10 6 is obtained from the integrated luminosity times the cross section of r-pair production, 



which is calculated in the updated version of KKMC [19j to be a TT = 0.919 ± 0.003 nb. 
We set the following 90% C.L. upper limits on the branching fraction products: B{t~ — > 
e"/o(980)) x B(/o(980) tt+tt") < 3.2 x 10~ 8 and B(r~ -> ^"/ (980)) x B(/ o (980) -> 

7T+7T-) < 3.4 X 10" 8 . 
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FIG. 4: Scatter-plots in the My - AE 1 plane with the ±20<7 area for (a) r~ — > e~/o(980) and (b) 
r~ — > /i~/o(980). The data are indicated by the solid circles. The filled boxes show the MC signal 
distribution with arbitrary normalization. The elliptical signal regions shown by a solid curve are 
used for evaluating the signal yield. No data events are observed in the elliptical signal regions. 



SUMMARY 

We have searched for lepton-flavor- violating r decays into a lepton and /o(980) meson 
using 671 fb _1 of data. No events are observed and we set the following 90% CL. upper limits 
on the branching fraction products: B(t~ -> e~/ (980)) x S(/ o (980) -> tt+tt-) < 3.2 x 1CT 8 
and B(t~ -> ^-/ (980)) x fi(/ (980) -> ir + n-) < 3.4 x 1(T 8 . This is the first search for 
these r LFV decay modes. 
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